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Abstract 

A simple collision model is employed to introduce elementary concepts 
of open system dynamics of quantum systems. In particular, within the 
framework of collision models we introduce the quantum analogue of ther- 
malization process called quantum homogenization and simulate quantum 
decoherence processes. These dynamics are driven by partial swaps and 
controlled unitary collisions, respectively. We show that collision models 
can be used to prepare multipartite entangled states. Partial swap dynam- 
ics generates VK-type of entanglement saturating the CKW inequalities, 
whereas the decoherence collision models creates GH Z-type of entangled 
states. The considered evolution of a system in a sequence of collisions is 
described by a discrete semigroup £1, . . . , £„. Interpolating this discrete 
points within the set of quantum channels we derive for both processes the 
corresponding Linblad master equations. In particular, we argue that col- 
lision models can be used as simulators of arbitrary Markovian dynamics, 
however, the inverse is not true. 



1 Open system dynamics 

The goal of these lectures is to present the elementary ideas and features of 
dynamics of open quantum systems by analyzing the properties of the simplest 
collision model we can think of. The material is based on papers [T] [31 131 13] on 
collision models coauthored by authors. It is not considered as a review paper 
on collision models or open system dynamics. 

It turns out it is surprisingly useful in physics to distinguish the concepts of 
isolated and open systems. The concept of isolated systems represents a simpli- 
fication of physical reality, which serves as a playground for a clear formulation 
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of elementary physical principles. It is postulated that the dynamics of isolated 
quantum systems is driven by Schrddinger equation [5] 

*^^^t = ^ti^t , (1) 

where Ht is a hermitian operator called Hamiltonian. As a consequence it 
follows that the state transformations are unitary, i.e. ~^ V'*' = Ut^t'tpt, 
where Ut^t'Ul^^, = U^^^,Ut^t' = I- If Ht is time- independent, i.e. Ht = H, 
then Ut^t' = exp [— — i)]. We can define Ut = exp{—jrHt) and write 

Ut^t' = Ut'Ul i^t = Ut^o- 

Interactions between a system under consideration and its environment are 
responsible for the violation of system's (dynamical) isolation. The joint evolu- 
tion of the system and the environment will be still unitary, however, the system 
itself undergoes a different dynamics. The point is that these interactions, be- 
cause of the complexity of the environment, are typically out of our control. 
And our wish is to invent models of system-environemt interactions capturing 
faithfully the key properties of the dynamics of the system alone. For a general 
discussion on the models of open system dynamics we refer to O [71 18] . In what 
follows we will restrict our analysis to a particular toy model of the open system 
dynamics. 

In the rest of this section we will introduce the model and elementary prop- 
erties of quantum channels. In the following two sections we will discuss ther- 
malization and decoherence processes within this model. In Section IV we will 
focus on entanglement created in the discussed collision models and finally, in 
the Section V we will derive master equations for these collision processes. 

1.1 Simple collision model 

Let H be the Hilbert space of the system of interest. States are identified with 
elements of the set of density operators (positive operators of unit trace) 

S{H) ={g:g>0, tr[g] = 1}. (2) 

We will assume that environment consists of huge number of particles, each of 
them described by the same Hilbert space T-l. In a sense, it forms a reservoir 
of particles. Moreover, we assume that initially they are all in the same state 
^, i.e. the initial state of the environment/reservoir is w = ^^JV each time 
step the system interacts with a single environment's particle. The interaction 
is described by a unitary operator U acting on joint Hilbert space 'H®'H. Thus, 
the one-step system's evolution is described by a map 

q' ^£[Q]^iv,^,[U{Q®OU^], (3) 

where g is the initial state of the system and trenv denotes the partial trace over 
the environment. In each time step a system undergoes a collision with a single 
particle from the reservoir. In a realistic scenarios these collisions are happening 
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Figure 1: A simple collision model. 

randomly, but we will assume that each reservoir's particle interacts with the 
system at most once (see Fig. [1]). Thus, the system evolution is driven by a 
sequence of independent collisions. After the nth collision we get 

QO-^ g'n^ trenv[C/(f?n-l ® = £[gn-l] = ^ "[po] • (4) 



1.2 Quantum channels 

The mapping S defined on the set of all operators /^("H) satisfies the following 
properties: 

1. £ is linear, i.e. £[X + cY] = £[X] + c£[Y] for all X,Y e C{n) and ah 
complex numbers c. 

2. £ is trace-preserving, i.e. tr[5[X]] = tr[X] for all operators X with finite 
trace. 

3. £ is completely positive, i.e. {£ > O for all positive operators 

Due to Stinespring dilation theorem the inverse is also true. If a mapping £ 
satisfies above three conditions, there is a unitary operator U and a state ^ such 
that the Eq. ([3]) holds. Such completely positive trace-preserving linear maps 
we call channels. 

For any state ^ G S{H) there exist a unit vector ^ T-L ® % such that 
^ = tr2|5')(^'|, where tr2 is the partial trace over the second subsystem. We call 
^E* a purification of ^. It is straightforward to verify that the following identity 
holds 

£[Q\^iT,n.[U{Q®OU^]=^^onA{U(^I){Q®\^)miU(^I)^]. (5) 

Let ipj is an orthonormal basis on % ® % and define linear operators Aj = 
{^pj\{U (g) 7)^) acting on the system. Then the right-hand side of the above 
equation reads 
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By definition the so-called Kraus operators Aj satisfy the normalization ~ 
/. Any mapping of the Kraus form £[g] = '^j Aj gAj for arbitrary set of op- 
erators {Ai,A2, . . . } fulfilling the normalization constraint determines a valid 
quantum channel. 

Channels with a single Kraus operator £[g] — AgA'' are unitary channels, 
because the normalization A^ A = I implies that A is unitary. Let us note that 
Kraus representation and also the number of Kraus operators are not unique. 
For more details we refer to [SI [TU] ■ 

2 Quantum homogenization as an analogue to 
thermalization 

The 0th law of thermodynamics postulates that an object brought in a contact 
with a reservoir of temperature T will equalize its temperature with the reser- 
voir's temperature. This process is known as the process of thermalization. Our 
attempt is to design a quantum analogue of this process. 

The concept of temperature is the key ingredient of the thermalization pro- 
cess. In fact, as a consequence it is the temperature that describes a state of 
the system, that is, it contains the relevant part of the information needed for 
the investigation of the system's properties. What is the quantum analogue of 
the temperature? There are several ways how this concept can be introduced. 
One can use the connection between temperature and the entropy, or energy per 
single particle. By choosing such analogue, the thermalization can be defined as 
a process that leads a system with a given temperature Ts to the temperature of 
the reservoir T^. In classical theory this process is achieved because of mutual 
interactions/collisions between the particles. In the process of collision they 
mutually (partly) exchange their microscopic properties (energy, momentum). 
Subsequently, the temperature of the system being in a contact with the reser- 
voir is changing and finally (approximately) equalize T^. Let us note that we 
have used the temperature as a property that can be attributed even to a single 
particle, however, such definition is meaningful only in the statistical sense. 

In quantum theory the state of a quantum system is described with a density 
operator. This operator contains all the information. Therefore, we choose 
single-particle density operator to play the role of the quantum temperature in 
our quantum analogue. So the goal is to design a model of the dynamics of the 
system-f reservoir such that the initial state of the system particle gs evolves 
into the state of the reservoir particles ^. Because of the property that finally 
all particles are described by the same state ^ we will refer to this process as 
to homogenization instead of the thermalization. In fact this notion is more 
appropriate to characterize the features of our model. 

Ideally the homogenization should implement the following transformation 

£•5 «) C ® •■■ ® C '-^ «) C «)•••«) C , (7) 

for all g and ^. For finite reservoirs this process will perform N ^ N + 1 
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cloning transformation, which would violate the quantum no-cloning theorem 
pn [T2] . Therefore, this process cannot be achieved in any dynamical model 
respecting the quantum theory. We will relax our requirements and investigate 
whether the homogenization can be implemented approximately via a sequence 
of (independent) collisions. Moreover, if initially g = meaning there is nothing 
to homogenize, then the evolution is trivial. Let us summarize our assumptions 
and problem. 

Definition 1. We call an interaction U a generator of J-homogenization S > 
is the following conditions are satisfied 

• Trivial homogenization. The joint state of the system and the environment 
after nth interaction reads 

n„ =[/„••• Uiigs <E> e^")C// • • • Ul , (8) 

where Uj = Us,j ^env\{j} describes the jth interaction (/env\{j} is the 
identity operator on the environment excluding the j th particle) . The for 
all C e S{H) 

^5 (g) ^ • • • (g) ^ ^ 17„ = ^5 (g) ^ (g) • • • ^ . (9) 

(n) 

• Convergence. Let gg be the state of the system after nth collision. Then 
there exist Ns such that for all n > Ng 

DigP,0<S. (10) 

• Stability of the reservoir. Let be the state of the jth reservoir's particle 
after the its interaction with the system. Then 

mpO<s, (11) 

for all j. 

Our task is simple: find U satisfying the above properties. 
2.1 Trivial homogenization 

Let us denote hy 'H± the symmetric and antisymmetric subspaces of "H g) 
A vector $ belongs to 'H± if 5$ = ±$, where S is the swap operator defined 
as S' = J2jk I'Pj ^ 'Pk){Vk <8) (pj\ for arbitrary orthonormal basis {tpj} of "H. It 
follows that for all g,^ e S{'H) 

S{g^OS^ ^^(^g- (12) 

Let us denote by P± projections onto symmetric and antisymmetric subspace. 
Then S' = P+-P_. 
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The condition of trivial homogenization implies that 



i7|V'0 V') = I'^i'fS) V') 



(13) 



for all unit vectors ip G H. This determines the action of the transformation 
U on the symmetric subspace H+. However, it docs not give any constraint 
on the antisymmetric subspace H-. Since % ®'H = H+ © H- we can write 
U = e^'^I+ ® V- , where V- is arbitrary unitary operator on the subspace H- and 
is the identity operator on the symmetric subspace Hj^. In fact, /+ = P+\'u+- 
Any of the coUisions of the form U = e^'^I+ © V- fulfills the conditions of trivial 
homogenization. 

Let us consider the case of d = dim?^ = 2 (qubit), for which the antisym- 
metric space is one dimensional. In one-dimensional Hilbert spaces the unitary 
operators are complex square roots of the unity, i.e. V- = e'^/_. Explicitly the 
transformation can be expressed as 



U\00) 

u\n) 

C/(|01) + |10)) 
C/(|01)-|10)) 



e'^lOO) , 
e'T|ll) , 

e*^(|01) + |10)), 
e»(7+/3)(|oi)_ |io)) 



Using the identities e}^ = e^('>'+'3/2)e-»/3/2^ = gi(j+P/2)^if3/2 discard- 

ing the irrelevant global phase factor e^'^''^^^^^ we get a one-parametric set of 
solutions 



= cos 77/ -|- i sin 775 , (14) 

where we set r] = —(3/2. 

In case of more dimensional space (where d > 2) the set of unitary trans- 
formations is larger. However, when we fix V+ = e'^/_, we get the same one- 
parametric set of solutions 

Ur, = COS T]I + i sin r]S = e'"^ . (15) 

Any element of this class of unitary operators we will call a partial swap. In 
what follows we will focus on the collision dynamics generated by partial swap 
coUisions. 



2.2 Pcirtial swap collisions 

In what follows we denote sin = s and cos r] = c. In the process of homogeniza- 
tion, the system particle interacts sequentially with one of the N particles of 
the reservoir through the transformation Un = cl + isS. After the first collision 
the system and the first reservoir particle undergo the evolution 

4"' ® ^ ^ Uriigf ® = c^gf ® C + ® + ics[S, gf <E> ^] . 
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The states of the system particle and of the reservoir particle are obtained as 
partial traces. Specifically, after the system particle is in the state described by 
the density operator 

9i'^-c'9f+s'^ + ^c.s[^,9f], (16) 
while the first reservoir particle is now in the state 

e, = s'gf+c'^ + ^cs[gf,^]. (17) 

When we recursively apply the partial-swap transformation then after the in- 
teraction with the rt-th reservoir particle, we obtain 

4")=c2 4"-i)-f.2^ + ^c.K,4"-i)], (18) 

as the expression for the density operator of the system particle, while the n-th 
reservoir particle is in the state 

2.3 System's convergence 

Our next task is to analyze the remaining conditions on the homogenization 
process. First task is to show that gi^"'' monotonically converges to a 5-vicinity 
of ^ for all parameters rj 0. We will utilize the Banach fixed point theorem 
|13j claiming that iterations of a strictly contractive mapping converge to the 
single-point contraction. This point is unique and determined by the condition 
£[^] = 5. A channel £ is called strictly contractive if it fulfills the inequality 
£)(£:[£.], f[C]) < kD{g,C) with < fc < 1 for aU £», ^ e S. 

First of all, it is easy to verify that ^ is a fixed point of £^ , where is defined 
by the Eq. (fT6|) . As a distance we will use the Hilbert-Schmidt distance 

D{gi,g2) = \\gi-g2\\-- tr[gl] + trigj] - 2tr[^)l^)2] (20) 

where gi, Q2 represent quantum states. Consequently, 

D^{£[g^],£[g2]) = c^D^g^, g^) + c^s^^, g^ - g2]\\'' (21) 

where we used that £[g\ = c^g + s'^£,+ics[£,, g]. Defining an operator A = gi — g2 
our aim is to prove the relation ||[^, A]|p < ||A|p. Assuming the spectral 
decomposition of ^ equals ^ = J2k ■^k\k){k\ and performing the trace in the 
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basis \k) we obtain 

||^A-A^||2 = tr[(^A-AO^(eA-Ae)]=2tr[CW]-2tr[(CA)2] 

= 2l]/'01A'li) - 2^. ^A^Afe|01A|fc)|2 
= . ,A|01A|fc)(fc|A|i) - 2^ ._^A,A,|01A|fc)|2 

= E j2A,2-2A,Afc)Ki|A|fc)p 

= E,,,(^' + 2A,A,)|0-|A|fc)|2 = ^(A, - Xkflijimi' 

< ^ \{j\A\k)\'=tv[A^A] = \\AW' 

where wc employed the hermiticity of A (|(j|A|fc)p = |(Z|A|7)p) and positivity 
of ^ which means |Ai — A2I < 1. The proved inequahty enables us to write 

D{£[e^],£[Q2])<\c\D{Qi,e2) (22) 

i.e. the contractivity coefficient k is determined by the parameter c — cos rj of 
the partial swap and the map £ is contractive whenever |c| < 1. This result is 
important because it ensures that in the limit of infinite steps the system will 
be not only (5-homogenized, but will be described exactly by the state ^ of the 
reservoir. It means the distance D{Qg'\^) vanishes with the number of 
interactions n. 



2.4 Stability of the reservoir 

In order to satisfy the last condition (stability of the reservoir) we need to 
evaluate the distances D{£^'p£). For each value of 77 we can specify 5 and Ns- 
However, the question is whether for arbitrary value of (5 > there is a suitable 
collision Ur, and how large is the value of Ng for which D(g^g'''\s,) < S. 

For the system we can combine the contractivity bound £>(f [^], f [^]) < 
\c\D{g,^) with the fact that £[^] = ^ to obtain an estimate 

D{^, gi"^) < \crD{^, qP) = V2\cr < 6 . (23) 

Consequently, 

but the potential values of 5 we need to specify from the stability of the reservoir. 
A direct calculation gives the following bound 

D{^,0 = u-c'^-s'e^r"^-ics[gf-'\m 

< s'\\^-gt-'\ + \cs\.\\[g^r"\m 

< s'D{^,g^r"^)+2\cs\-urH-m\ 

< y2|sc|(|s| • |cr-2 + \/2) 

< V2\sc\{l + V2) = 5 . (25) 
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Important is that S can be arbitrarily small, hence (5-homogenization is achiev- 
able for a restricted class of collisions Ur, satisfying the identity (2-|-\/2) | sin r] cos r/l < 
6. The number of steps Ns we achieve if the value of S is inserted into the formula 
for Ns, i.e. 

^ ln[(l + 72)|sc|] 

" hTM ■ ^ ^ 

Let us note that one can derive more precise expression for Ng and S, but this 
is not important for our further purposes. 

In summary, the class of partial swap operators Ur, satisfies the conditions 

for (5-homogenization process. 

Remark 1. Swap in the classical theory. Let us note that in the classical picture 
of thermalization process we also use model of mutual collisions of the system 
with its thermal reservoir. In these (inelastic) collisions the particles exchange 
energy and momentum. In some sense, this process can be understood as a 
transformation that partially "swaps" energy and momentum of involved par- 
ticles. Thus, from this perspective the homogenizing properties of partial swap 
collisions are not that surprising and can be viewed as the quantum analogue 
of classical collisions. 



2.5 Invairiance of single-peirticle average state 

In this section we will allow particles forming the reservoir interacts among each 
other, but interactions are always pairwise. Partial swap operators possess one 
unique property, which is important from the point of view of more realistic 
thermalization process and reservoir's stability. Consider O being an n-partite 
composite system. It can be written in the form 

O = ® • • • ® e„ + r , (27) 

where F is a traceless hermitian operators and qj = tTjQ are the states of 
individual subsystems. Applying a partial swap collision between jth and kth 
subsystems we obtain a state 

a' = c^n + s'^SjknSjk + ics[n,Sjk] (28) 
= e'i0---0^; + r', (29) 

where Sjk denotes the swap operator between the subsystems j and k, 

g'j = c^Qj + s'^gk+ics[gj,gk\ (30) 
Qk = c^Qk + s'^Qj - ics[gj,gk] (31) 
g'l = gi for l^j,^, (32) 

and r' is the traceless part. It is straightforward to see that 

^one = ~ H = ~ = ^one , (33) 

j j 



9 



thus, in the cohision model provided by partial swap collisions the average one- 
particle state 'Qanc preserved. If taken the single-particle state as the quantum 
analogue of temperature, then this property means that partial swap interac- 
tions preserves the "quantum" temperature. 

3 Quantum decoherence via collisions 

In this section we will focus on quantum dynamics of decoherence. Let us note 
that in literature on quantum information the decoherence is used to name any 
nonunitary channel and the channels we are interested in are called (generalized) 
phase-damping channels. In accordance with our previous work [3| we will 
understand by decoherence a channel with the following properties: 

• Preservation of the diagonal elements of a density matrix with respect 
to a given (decoherence) basis. Let us denote hy B = {ipi, . . . , ipd} the 
decoherence basis of H. Then 

£ o diagg = diagg , (34) 

where diagg is a channel diagonalizing density operators in the basis B, 
i.e. g — diagg [g]. 

• Vanishing of the off-diagonal elements of a density matrix with respect to 
a given (decoherence) basis. In particular, asymptotically the iterations 
of the decoherence channels results in the diagonalization of the density 
matrix, i.e. 

lim f " = diagg . (35) 

n— >oo 

Our tasks is similar as in the previous part; find collisions U generating the 
decoherence in basis B. 

The basis preservation implies that 

1^, ^e*"^'"!^, ®(^;.fc), (36) 

where rjjk are phases and for each j the unit vectors {</3jj.}fc form an orthonormal 
basis oiH. It is not difficult to verify that the unitary operators describing the 
collisions have the following form 

C/ = ^|^,)(^,|®V,, (37) 
i 

where Vj — J2k^^^^'°\'fi'jk)iVk\ are unitary operators on individual systems of 
the environment. Such operators form a class of controlled unitaries with the 
system playing the role of the control system and the environment playing the 
role of the target system. 
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A single controlled unitary collision induce the following channel on the 

system 

£[g\ = trenv[C/(e® Ot^l = E l^Msl^kHVkMVj^V^] . (38) 

Let us note an interesting feature that under the action of this channel the 
matrix elements of density operators do not mix, i.e. the output value of Qj^. = 
{ipj\£[g\ipk) depends only on the value of Qjk = {ipj\gipk)- In order to fulfill the 
second decoherence condition it is sufficient to show that for j ^ k 

\Q'jk\<\Qjk\- (39) 
For any unitary operator W = J2i e'"''|Z)(/| 

I I 

where pi = are probabilities. This inequality is saturated only if ^ is an 

eigenvector ofW. It follows that the inequality 

141 = l^ifcl • MVj^Vk]\ < \0jk\ (41) 

holds, hence the off-diagonal elements are non-increasing. They are strictly 
decreasing if and only if |tr[V^-^y^^]| < 1 for all j ^ k, which is the case if ^ is 

not an eigenvector of any of the operators V^Vj. 

In summary, we have shown that decoherence processes are intimately re- 
lated with controlled unitary operators. In particular, the decoherence processes 
are generated by controlled unitary collisions U = \ ® Vj satisfying 

|tr[yj-^yj]| < 1 for all j ^ k. In such case 

lim £^[q] = diag^l^] , (42) 

n— >oo 

where diagg is a channel diagonalizing the input density operator in the basis B. 
Let us note that there arc always states ^ (eigenvectors of V^Vj) for which the 
decoherence is not achieved, because for them some of the off-diagonal elements 
are preserved. On the other side, for each controlled unitary collision U there are 
always states of the environment ^ inducing the decoherence of all off-diagonal 
matrix elements. 

3.1 Simultaneous decoherence of the system and the en- 
vironment 

Could it happen that both the environment and the system arc decohering 
simultaneously? In the same decoherence basis? We will see that answers to 
both these questions are positive. 
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For the considered collision model driven by controlled unitary interactions 
we obtain that the evolution of the particles in the reservoir is described by the 
following channels 

AA[e] = tvsysiUig ® OU^] = J2 QnVj^Vj . (43) 

j 

Such channels are known as random unitary channels. 

It is known that random unitary channels can describe decoherence channels 
and it was a surprising result that decoherences are not necessarily random 
unitary channels [15]. We do not need to go into the details of this relation to 
see that both the system and the environment can decohere. It follows from our 
previous discussion that environment will decohere if U is controlled unitary 
transformation of the form 

U^J2U,^\^,}{^,\, (44) 
i 

where {tpj} is the decoherence basis of the particles in the environment. That 
is, the question is whether there are unitary operators that can be written as 

U = J2U,^ IV-,) I = 1^^) 1 ^ ' (45) 

j 3 

for suitable decoherence bases {(/Sj}, {"0^} and unitary operators Uj,Vj. 
Example 1. CTRL-NOT. 

f/ctrl-NOT = |0)(0|®/+|l)(l|®a, =/®|+)(+|+f7,®|-)(-| (46) 
= |0+)(0+| + |0-)(0-| + |l+)(l + |-|l-)(l-|, (47) 

where |±) = (|0) ± |l))/\/2. 

Proposition 1. A unitary operator U describes a collision with simultaneous 
decoherence of both interacting systems if and only if U = \ipj){ipj \ ® Vj and 
the unitary operators Vj commute with each other. 

Proof. By definition, the simultaneous decoherence requires simultaneous preser- 
vation of decoherence bases {(pj} and {^k}i i.e. 

\ipj (g) Vfe) ^ e'"^^" (E) i^k) ■ (48) 

It follows that 

C/ = ^e^''^-|¥'.,-)(^.l®l0O(V'fc|, (49) 

jk 

and we can define 

Uk^Y. I ' = E e^'^" IV'fe) i^k I , (50) 

3 k 
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to get the required expression 

U^Y.Uk^ I Vfc) (^fe I = ^ \^Vj. (51) 
k j 

It is straightforward to see that operators {Vj} and {Uk} form sets of mutually 
commuting elements. 

If the unitary operators Vj commute with each other, then in their spectral 
form V, = e'"^^ IV'fe) (^fe I- Clearly, U = i^j I ® = Efe Uk ® | Vfe) (^fe I 

with Uk = J2j Therefore, U generates decoherence in both inter- 

acting systems. □ 

Setting ipk = fk the system's decoherence basis and the environment's de- 
coherence basis coincide. In such case the system and the environment will 
decohere with respect to the same basis. An example of such interaction for the 
case of qubit is the CTRL-Z operator 

C/ctri-z = |0)(0| ® / + |1)(1| ® a, = / ® |0)(0| + (T, ® |1)(1| , (52) 

generating decoherence in the computational basis |0), |1). Moreover, this col- 
lision is symmetric with respect to exchange of the role of the control and the 
target. 



4 Entanglement in collision models 

The phenomenon of quantum entanglement is considered as one of the key prop- 
erties of quantum theory. Our aim is not to discuss the philosophical background 
of this concept [16| I17j , but rather focus on the dynamics of entanglement in 
the considered collision models. We say a state is separable if it can be written 
as a convex combination of factorized states, i.e. 

ftcp-5]p,^P^$5---^C«. (53) 

j 

A state is entangled if it is not separable. The question whether a given state 
is entangled, or separable turns out to be very difficult. Partial results are 
achieved for the system of qubits and therefore in this section we will restrict 
our discussion only to qubits. However, some of the qualitative properties will 
be extendible also to larger systems. 

A tangle is a measure of bipartite entanglement defined as 

^(^) = ^ min Vpjr(V'j) , (54) 

where t(V') = 5ii„(tri|V')(V'|) = 2(1 - tr[(tri|V')(V'|)']) = 4dettri[|V')(^|] is 
the linear entropy of the state of one of the subsystems. Let us note that for 
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pure bipartite states 5'iin(tri|^)(V'|) = 5'iin(tr2|V')('!/'l)- Due to seminal paper of 
Wootters \18\ t{uj) — C^(w), where C(uj) is the so-called concurrence 

C = max{0, 2 maxlV^} - ^ ^A"} , (55) 

j 

where {Xj} are the eigenvalues of i? = ^{<^y ^ ay)uj*((Ty ® ay). 

For multiqubit systems the tangle satisfies the so-called monogamy relation 
(originally conjectures by Coffman et al. 19 and proved by Osborne et al.|20p 

where Tj = t^j and j denotes the set of all qubits except the jth one, i.e. Tj is a 
tangle between the jth qubit and rest of the qubits considered as a single system. 
The above inequality is called CKW inequality. If the multiqubit system is in a 
pure state 5" then 

Tj- =4dettrj[|*)(1'|]. (57) 

Example 2. Saturation of CKW inequalities. At the end of the original paper 
jl9j it is stated that states in the subspace covered by the basis {|l)j<8i |0®^~^)j} 
saturates the CKW inequalities. In what follows we will show that any state of 
the form 

N 

\^) = aolon +E"^ |1)^- ® |0^(^-^)) (58) 

saturates the CKW inequalities. 

The reduced bipartite density matrices are of the form 



(59) 



and only one element determines the tangle of such state, namely 

t{q) - 4//* . (60) 

It follows that only the matrix element standing with the |01)(10| term will be 
important for us. Direct calculations lead us to value / = OLjO*). for the pair of 
jth and fcth, hence 

r,fe =4|a,f lafcp. (61) 

In the next step we evaluate the tangle between the j-th particle and the rest 
of the system. The state of single particle is described by matrix 
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d* 
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e* 
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Now it is easy to check that 



r, = 4 det = |a,f ^ |2 ^ ^ r,, (63) 

and therefore the CKW inequahties are saturated, i.e. Aj = tj — J2k^j '^ifc — ^ 
for all j. 

4.1 Entanglement in partial swap collision model 

A general pure input state equals 

\'i/)^\ij)s<E)\<f)i(E)---(E)\ip)N- (64) 

Since [/, ^ (8) T^] = [5, F (8) 1^] =0 for all unitary operators F on "H it follows 
that also [C/^, FiSiI^] =0 for all unitary V. Moreover, the value of entanglement 
is invariant under local unitary transformations. Therefore, without loss of 
generality we can set \(p) — |0) and jV') = a|0) + After n interactions 



= a|0®(^+i)) + /3c"|l)5 |0®^) + |l)z |0®^T) 



1=1 



where Nj denotes a system of N qubits obtained by replacing the Zth qubit from 
the reservoir by the system qubit. 

Since |5'„) belongs to the class of states discussed in Example [21 we know 
that the state saturates CKW inequalities. We can used the derived formulas 
to write 

r.kin) = 4|/3|4s4c2(^-+'=-2) ; (65) 
rokin) = 4|/3|S2c2("+'=-i) ; (66) 
T,{n) - 4|/3|Vc2(j"-i)(l-sV(J-i)); (67) 
To(n) = 4|/3|V"(l-c^"). (68) 

Let us note that these formulas are valid only if j, k < n. Otherwise the quan- 
tities vanish. 

These results show that the system particle acts as a mediator entangling the 
reservoir particles which have never interacted directly. It is obvious that the 
later the two reservoir qubits interact with the system particle, the smaller the 
degree of their mutual entanglement is. Nevertheless, this value remains con- 
stant and does not depend on the subsequent evolution of the collision model 
(i.e., it does not depend on the number of interactions n). On the other hand 
the entanglement between the system particle and kth reservoir particle (k arbi- 
trary) increases at the moment of interaction and then monotonously decreases 
with the number of interaction steps. 
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Figure 2: Entanglement in the homogenization process. 



4.2 Entanglement in controlled unitary collision model 

As before, let us assume that initially the system and the reservoir of N qubits 
are described by a pure input state 

\n) = \i;)s 1(E) ■■■(E) \ip)N , (69) 

with ip — a\0) + /3\1). After n controlled unitary collisions we get 

la.) = [a|0) E) \(pT') + m E) {(fT)] ^ |^®(^-")) , (70) 

where we set \(po) — Vq\ip) and \Lpi) = Vi\(p). 

Tracing out the system we end up with the reservoir described by the state 

.._(n) = {\a\'\<pr){^r\ + wi'iffnifri) ® , (71) 

which is clearly separable, i.e. Tjk{n) = 0. In this case no entanglement is 
created within the environment. 

A state of the system qubit and fcth qubit from the reservoir equals (provid- 
ing that n > k) 

Qokin) = |anO(^o)(0<^o| + |/3ni'/'i)(l'^i| 

+a/3* I (^0 l^i) 1^"-') |0^o> {l(Pi I + c.c , (72) 

thus for the tangle we have 

ro,(n) = i\af\/3f\{Vo\^i)f^''-'^\{VoU)\' ■ (73) 
If n < k then TQk{n) — 0. Further, 

go{n) = |ano)(0| + |/3ni)(l| + ar|(^ol^i)nO)(l| + c.c., (74) 
Qkin) = |an^o)(^o| + |/3|'|^i)(^i|, (75) 

(76) 

resulting in the following values for tangle 

ro(n) = 4|an/?|2(l-|(^o|^i)|'"), (77) 
Mn) = 4\am'\{^oU)\\ (78) 



16 



1 

. 8 

. 6 
. 4 
. 2 



2 4 6 8 10 12 14 n 

Figure 3: Entanglement in the decoherence process. 



In summary, during the decoherence collision model the qubits in the reser- 
voir are not entangled at all. Interaction of the system qubit with an individual 
reservoir's qubit entangles these pair of qubits. However, with the number of 
collision this bipartite entanglement vanishes and finally Tok — for all reser- 
voir's qubits k. The entanglement between a fixed qubit of the reservoir and rest 
of the qubits preserves its value it gathered after the interaction. The entangle- 
ment between a system qubit and the reservoir is increasing with the number 
of collisions and approaches the maximal value to(oo) = 4|ap|/3p. 

5 Master equations for collision models 

Within the collision model the system undergoes a discrete dynamics 

g ^ £i[g] ^ £2[g] ^ > £n[g] ^ ■ ■ ■ , (79) 

where £k — £ o ■ ■ ■ o £ = £''. Let us define £o as the identity map I. It 
is straightforward to see that the sequence of channels £q,£i,£2, ■ ■ ■ form a 
discrete semigroup satisfying the relation 

£n ° £m = £n+m (80) 

for all n, TO = 0, 1, 2, ... . In this section we will investigate the continuous inter- 
polations of such discrete semigroups. Interpreting the continuous parameter as 
time we will derive first order differential master equation generating the con- 
tinuous sequence of linear maps. Our derivation of the master equation will be 
purely heuristic. We will simply replace n in £„ by a continuous time parameter 
t and after that we will have a look what are the consequences. 

Let gt = ft[po]- The time derivative of the state dynamics results in the 
following differential equation 

where Qt is the generator of the dynamics £t. Let us note that this approach 
is very formal. There is no guarantee that the mapping £t is a valid quantum 
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channel for each value of t and that the set £t is indeed continuous. Another 
problematic issue is the assumed existence of the inverse of the map for ah 
t>0. All these properties need to be checked before we give some credit to the 
derived master equation. 

5.1 One-parametric semigroups 

In a special case when Qt is time-independent the solution of the differential 
equation 

with the initial condition gt=o — go can be written in the form 

gt^e^*[go]; £t = e^K (83) 

Since e^(*+') = e^"e^* it follows that 

£to£s= £t+s (84) 

for all t, s > and £q = I. In such case the continuous set of linear maps 
{£(} form a one-parametric semigroup, which are playing important role in the 
theory of open system dynamics. Semigroups of channels are often used to 
approximate real open system's evolutions. In particular, they are satisfying 
the so-called Markovianity condition^^ restricting the memory effects of the 
environment. The following theorem specifies the properties of the generator Q 
assuring the validity of the quantum channel constraints on £t — e*^ . 

Theorem 1. (Lindhladj2^) £t — e^* are valid quantum channels for all t > 
if and only if Q can he written in the form 

Q[Q\ = -\[H,g] + \Y. c,fc([A,e, Afc] + [A,, ^,Afe]) , (85) 

jk 

where {/, Ai, . . . , Ajj2„]^} is an orthonormal operator basis consisting of hermi- 
tian traceless operators Aj, i.e. tr[AjAfc] = 6jk, the coefficients Cjk form a 
positive matrix and H is a hermitian traceless operator on T-L. 



5.2 Divisibility of channels 

Any collision U induced some channel £. Our interpolation should smoothly 
connect this channel with the identity map X such that the connecting line is 
inside the set of channels. Since the set of channels is convex any two points are 
connected by a line containing only proper quantum channels, i.e. interpolation 
is surely possible although its continuity is still not guaranteed. Let us have a 
look on the possibility to interpolate X and £ such that £ = e'^^ for some time 
T and time-independent Lindblad generator. Without loss of generality we can 
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set T = 1, i.e. E — . The semigroup property implies that E — Ee o Ei-e for 
any < e < 1, or even stronger E = Ee^'^ ■ 

It is a surprising fact [23, that there are channels E, which cannot be expressed 
as a nontrivial composition of some other channels £i , £2 • Consequently not all 
channels can be part of some one-parametric semigroup. Before we will show 
one example, let us formally define the indivisibility. 

Definition 2. We say a channel E is indivisible HE — E10E2 implies that either 
El, or E2 are unitary channels (such decomposition is called trivial). 

Let us note that one of seemingly counter-intuitive consequences is that 
unitary channels are indivisible. In a suitable representation the composition of 
channels is just a product of matrices. Thus for determinants 

det {El o £2 ) = (det f 1 ) (det E2 ) • (86) 

Let us denote by Emm a channel minimizing the value of determinant, i.e. Emm = 
argmiuf detE. Assume Emin — Ei o E2- Since the value of detfmin is negative 
and for channels | detf | < 1, it follows that either detfi = 1, det £2 = det Emin, 
or det £2 — l,det£i — det But det£ — 1 implies f is a unitary channel, 
which proves the triviality of any decomposition of the minimum determinant 
channel. If d = 2 (qubit), then one of the minimum determinant channels is the 
optimal universal NOT^IS\. i.e. 

fmin[^^] = J(/+e^), (87) 

for which detfmin = —1/27. The mapping g —?' known as the universal NOT 
is not a proper quantum channel, because it is not completely positive and, thus, 
unphysical. Let us note that channels of minimal determinant are not the only 
indivisible channels. As far as the authors know a complete characterization of 
indivisible channels is not known. 

For our purposes the introduced concept of indivisibility is not sufficient as 
we are interested in the existence of continuous semigroup between I and E 
generated by a unitary collision U. We need E to be infinitely divisible, i.e. 
E = f jy^ for all n > 0. It was shown by Denisov[24] that infinitely divisible 
channel can be expressed in the form E = J- o , where J- is an idempotent 
channel (J-"^ = T) satisfying J^Q = TQT and is a Lindblad generator. We are 
interested in cases when T =T, thus, E = . Then it follows that the sequence 
of channels £" = e"^ determined by the collision model belongs to the one- 
parametric semigroup e*^. In other words a discrete semigroup is interpolated 
by a continuous one. We will see that this is the case for the collisions considered 
in the homogenization and decoherence processes. 

For a given channel induced by a single collision a question of interest is 
whether E = for some Lindblad generator Q. In other words whether Q = 
logE is a valid generator of complete positive dynamics. Since channels E are 
not necessarily diagonalizable, the notion of the matrix/operator logarithm is 
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not completely trivial. The logarithm is not unique and any operator X such 
that = f is called \og£. We need to search all logarithms in order to find 
a valid generator of the Lindblad form. Evaluating the logarithm will give us 
a time-independent generator (if it exists) interpolating the discrete semigroup 
= For both considered collision models the evaluation of the logarithm is 
not difficult. However, for illustrative purposes in our analysis of homogenization 
we will follow the ad hoc procedure exploiting (jSip , which (if successful) could 
capture also time-dependent generators. For the case of decoherence we will 
evaluate the logarithm of the induced channel. 

5.3 Bloch sphere parametrization 

The set of operators form a complex linear space endowed with a Hilbert- Schmidt 
scalar product {X, Y) = tr[X^y]. Let (/?!,..., (^^ be an orthonormal basis of H. 
Then the set of d^ operators {Ejk — \ipj){ipk\} form an orthonormal operator 
basis, i.e. ti[E'l^,^,Ejk] = {(pj'\(pj){ipk\'Pk') = ^jj'hk'- 

In what follows we will restrict to the case of two-dimensional Hilbert spaces 
% (qubits). This will be perfectly sufficient for our later purposes, because 
we will focus on qubit collision models of homogenization and decoherence. 
However, many of the properties and procedures can be easily extended to 
more dimensional case. 

Let us start with a very convenient operator basis of qubit systems - the set 
of Pauli operators 

/=|0)(0| + |1)(1|, a. = |0)(l|-|-|l)(0|, 
ay = -tm{l\-\l}{0\), a. = |0)(0|-|l)(l|, 

where {|0), |1)} is an orthonormal basis of 7i. These operators are hermitian, 
aj = crj, and mutually orthogonal, i.e. tr [ajCTfc] = 26 jt (we set ao = I)- They 
are also unitary, but this property will not be very important for us. 
Using Pauli operators the density operators takes the form 

g^^{I + f-a), (88) 

where a = {ax, cry,az) and f = tT[ga] is the so-called Bloch vector. The positiv- 
ity constraint restricts its length \r\ < 1, i.e. in the Bloch vector parametrization 
the qubit 's states form a unit sphere called Bloch sphere. Let us note that such 
simple picture of state space does not hold for more dimensional quantum sys- 
tems. 

A quantum channel £ is acting on Pauli operators as follows 

g' ^S[g]^ + x£[ax] + yS[ay] + z£[a,]) . (89) 

Using £[crj\ = J^k^kjO'k with £kj = ■^tr[o'i;£[crj]]. The tracepreservity ensures 
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that Sqo = ^tT[S[I]] = 1 and £oj = tr[£[aj]] = for j = x, y, z. That is, 



/ + I f fco + ^^^^1 I 

k—x,y,x \ j—x^y^x 



(90) 



Comparing with the expression = i (/ + r* • a) we get that under the action 
of a channel £ the Bloch vectors are transformed by an afSne transformation 



r r = t + Tf, 



(91) 



where tj = £jo and T is a 3x3 matrix with entries Tjk = £jk for j, k = x,y, z. 
The Lindblad generator takes the form of 4 x 4 matrix 



f 







g = 



\ 



310 .911 .912 313 

320 521 .922 323 
V 330 531 532 333 / 



(92) 



where gjk = ■^tr[ajg[ak]]- Inserting the Lindblad operator-sum form 



into the matrix expression we obtain 



9jk 

9k0 



I 



(93) 

(94) 
(95) 



The inverse relations express the parameters Cjk and hj via the elements of the 
matrix Q 



_ 932-923 L _ 913-931 l,^ _ 321-312 

"l — A , "•2 — / , "3 A 



- 9n ~ kY.k 9kk 



(96) 



Cl2 = 5(512 + 321 - «53o) , C21 = I (512 + 521 + «530) , 
C23 = 5(523 + 332 - »5lo) , C32 = 5(523 + 532 + »5io) , 
Cl3 = 5(513 + 331 + i.920) , C31 = 5(513 + 531 - 1520) • 

5.4 Master equation for homogenization collision model 

The collision model of quantum homogenization driven by partial swaps results 
in the sequence of channels being powers of a channel £ defined in Eq. (ll6|) . i.e. 



= (?r + s^w — csw x r , 



(97) 
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where w is the Bloch vector associated with the state f and c = cos 77, s = sin 77. 
Choosing a new operator basis Sj = VajV^ such that ~ ^{I + wS^) we get 





-CSW 









csw 

^2 







(98) 



Let us note that the basis transformation aj — >■ Sj 



corresponds to a rotation of 
the coordinate system of the Bloch sphere representation. Moreover, we used 
that the partial swap commutes with unitaries of the form V (E) V, hence the 
form of the induced map ([W)) is unaffected, only the vectors are expressed with 
respect to different coordinate system. 

Let us introduce an angle 6 = arctan(z«s/c) — arctan(u' tanr/) and a pa- 

1/ ' ' 



rameter q 



■w'^s'^. Using the identity cos arctan(a;) 



qcos6 = cos 77 = c and gsin6' = w sin 77, thus 



/ 


1 








' \ 







cq cos 9 


cq sin 9 










^cq sin0 


cq cos 9 





V 


s w 











'1 + x'' we get 



(99) 



The reason for such parametrization becomes clear if we evaluate powers of £ 



( 







c^q^ cos n9 
—c^q^ sin n9 



q 









sr 
cosrj6 




In the next step we make ad hoc assumption and replace n by t/r, where r 
is some time scale and i > is a continuous time parameter. Introducing the 
parameters 

we end up with the continuous set of trace-preserving linear maps 



(100) 







.-rot 



-r,t. 



cos fit 



-Tit 









e"^'"* sin nt 
e~^2*cosf7i 








-rit 



(101) 



It is not difficult to see that they form a one-parametric semigroup {£t o £g = 
£t+s)^ but one needs to verify whether for each t these maps are completely 
positive. Therefore, we calculate the generator Q and verify its properties. In 
particular, we obtain a time- independent generator 



/ 











' \ 







-La 


-n 










n 







v 
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Using the Eas. ([M| we can verify that the parameters /ii = = 0, = — 
are real and that the matrix 



-Ti/2 -i2wTi 
C= I i2wTi -ri/2 I (102) 

-(ri + 2r2)/2 

is positive, hence, the generator Q is of Lindblad form. In the operator-sum 
(Kraus) form the master equation reads 



dg . D, 

'dt ~ ^2h 



^T^^iS,, g] - iwTi{<S>.,y[g] + <^>y,[g]) (103) 



r r -j- 2r 

where $jfe[gi] — ^{[cjg, <Tk] + [(Jj, g(Jk])- After a little algebra we obtain 



dg . Q 

'dt ^ 

1„ „ „ „ „ s 1 



— = i — [Sz,g]~iwri{SxgSy-SygSx+igS^+iS^g) (104) 



-TiiS^gS^ + SygSy - 2g) + ■^{'^^2 - ^i){S:,gS:, - g) . 

Let us note that for the special choice of parameters the above master equa- 
tion coincide with the master equation used to model the spontaneous decay 
of a two- level atom. In particular, setting Fi = 2T2 = 27 = — - In cos 77 and ^ 
being a pure state {w = 1) we get 

d n 7 

— g = -i'l^[Sz,9] + -jp-S-gS+ - S-S+g- gS+S-) (105) 
where we used S± — {S^ ± iSy)/2. 



5.5 Master equation for decoherence collision model 

In this section we will repeat the same steps as in the previous one, but for 
collisions described by controlled unitary interactions. Without loss of generality 
we will assume that the decoherence basis coincide with the eigenvectors of Sz = 
Va^V^ operator. The channel induced by a collision U = |0)(0|(8)Po + |l)(l|'8)Vi 
reads 

/I \ 

„ A cos 09 A sin 09 

£ = „ , . , ^ ^ , (106) 
— Asm(p \ cos ip ^ ' 

\ 1 / 

where tr[ViVoC] = Ae'^. We can follow the same derivation of master equation 
as in the case of homogenization [4j. However, for the illustration purposes 
we will describe the second procedure based on evaluation of \og£. Unlike the 
case of homogenization the decoherence channel £ defines a hermitian matrix. 
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thus, the logarithm is pretty easy to calculate exploiting the simple functional 
calculus for hermitian operators. 

Let us start with the observation that 



A(cos(^/ + isin(^(7j,) = ei"^-fe*('^+2'='^)^« = ^in A/+»(^+2fex)a„ 



(107) 



where k is arbitrary natural number. For fc = the logarithm is called princi- 
pal. For our purposes it is sufficient to consider only this one, because we are 
restricted to angles inside the interval [0, 27r]. Since this matrix is the central 
part of the matrix £, it follows that 



g = log£ 



/ 

In A -ip 

(p In A 

\ / 



(108) 



Gt 



is the generator of the semigroup dynamics containing the channel £. Is f t = e 
a valid quantum channel for any t7 

Using the relations specified in Eas. ([M| we found that the only nonzero 
parameters are 



/i3 = (/3/2, C33 = -lnA/2. 
Since C33 > it follows that the matrix C is positive. Therefore 



(109) 



Gig] 



In A 
In A 



{[Sz, gSz] + [SzQ, Sz]) 



{SzqSz - q) ■ 



defines a correct Lindblad generator. Using H = ^ipSz and 7 = (21nA)/(^^ we 
obtain a well-known master equation 



(110) 



which is used to model the decoherence also for more dimensional systems. 



6 Conclusions 

In these lectures we introduced and investigated the simple collision model to 
capture relevant features and properties of quantum open system dynamics. 
We focused on two particular collision models determined by the choice of the 
unitary transformations describing the individual collisions: 

• Homogenization induced by the partial swap interactions — cos rjl + 
isinrjS. This process was motivated by the classical thermalization pro- 
cess (0th law of thermodynamics) . 
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• Decoherence induced by the controhed unitary interactions U = \(pj) 
Vj . This process describes the disappearance of the quantumness of quan- 
tum systems. 

It is known that creation of entanglement requires interactions between quan- 
tum systems. In our collision model we start from completely factorized state. 
It is an interesting question what type of multipartite entanglement is created 
via a well-defined sequence of bipartite collisions. We haven't provided a defi- 
nite answer to this problem, hut we illustrated that the creatcid entanglement 
of the particular collision models is not trivial. In particular, the partial swap 
interaction creates W-type of entanglement for which all the involved systems 
are pairwisely entangled although they did not interact directly. Moreover, dur- 
ing the whole evolution the CKW inequalities are saturated. The entanglement 
created in the decoherence collision model is of completely different quality. In 
this case the created entanglement is of GH Z-typc meaning that pairwisc en- 
tanglement in the reservoir vanishes, however the system is still entangled. The 
reservoir itself (the system is traced out) is in a separable state. The collision 
models can be understood as a (simple) preparation processes aiming to create 
multipartite entanglement. And it is of interest to understand what quality and 
quantity of entanglement it is capable of. 

The well-defined sequence of collisions and specific initial conditions imply 
that the discrete time evolution of the system is described by a discrete semi- 
group of natural powers of S. May this this discrete set of channels be inter- 
polated by a single one-parametric continuum of channels? Is this continuum a 
semigroup of channels, or not? Luckily, we have derived that the interpolating 
continuum for both cases of homogenization and decoherences. By deriving the 
corresponding master equations and testing the Markovianity of the generator 
we showed that these sets form one-parametric semigroups of channels. Let us 
note that a collision model simulation of given semigroups of channels is easy. 
Each semigroup £t can be discretized by introducing a parameter r and set 
£n = Snr- Since the choice of an interaction U inducing is not unique, also 
the collision model is not unique. However, the particular collision dynamics of 
the system will be the same. 

Let us note an interesting fact. Fix a unitary collision U. Even if ^, ^' induce 
a channels £, £' with valid generators Q, Q' , the convex combination A^-|-(l — A)^' 
does not have to be associated with a correct Lindblad generator. There are col- 
lision models that cannot be interpolated by Markovian dynamics. Surprisingly, 
collision models provides richer dynamics than Lindblad's master equations. To 
sort out example consider Markovian models developed in order to describe true 
open system dynamics can be efficiently simulated (approximated) in collision 
models. We believe that such toy collision models provide an interesting play- 
ground capturing all the conceptual features of quantum open system dynamics. 
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